Optimal strategy for a single-qubit gate and trade-off between opposite types of 

decoherence 
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We study reliable quantum information processing (QIP) under two different types of environment. 
First type is Markovian exponential decay, and the appropriate elementary strategy of protection 
of qubit is to apply fast gates. The second one is strongly non-Markovian and occurs solely during 
operations on the qubit. The best strategy is then to work with slow gates. If the two types are both 
present, one has to optimize the speed of gate. We show that such a trade-off is present for a single- 
qubit operation in a semiconductor quantum dot implementation of QIP, where recombination of 
exciton (qubit) is Markovian, while phonon dressing gives rise to the non-Markovian contribution. 



Already at the early stage of quantum information 
theory, the implementation of quantum computing was 
found to be extremely challenging due to decoherence 
processes 0. It was claimed, nevertheless, that a quan- 
tum computer could reliably perform quantum algo- 
rithms once the error per operation is below some thresh- 
old level which is approximately S ~ 10~^ 01 (for some 
caveats coming from memory of environment see Q ) • In 
any case, minimizing the error per gate is a reasonable 
strategy. Usually, one considers the decoherence time 
of the system and the time of performing operation 
(quantum gate) Tg, so that the error is given by the ratio 
S — Tg/Td 0. This suggests to search for systems with S 
as small as possible, and the obvious strategy to diminish 
the error is to apply fast gates. One then tacitly assumes 
that the process of decoherence is independent of running 
the gate, which formally means that it is a Markovian 
process, where the error, to first order, grows linearly in 
time. Gate speed-up may be achieved by selecting ma- 
terials to provide favorable spectrum characteristics Q 
or by applying techniques reducing unwanted transitions 
within the register space 6] as well as outside this space 
(leakage) 7]. 

The assumption of Markovian character of noise is by 
no means obvious and, indeed, it has recently been found 
that due to non-Markovian effects it may prove to be 
completely invalid. The notion of minimal decoherence 
model was introduced where the error occurs solely during 
gate operation, and it grows for fast gates a,sd ^ 1 /r^ (for 
the spectral density of the reservoir ~ w^). An example 
of such a model is a degenerate system interacting with a 
bosonic field at vacuum via dipole interaction. This kind 
of error favors the counterintuitive strategy of slowing 
down the gate operation. 

In real systems, this type of decoherence competes with 
the Markovian damping, although this effect may be cov- 
ered by other errors (e.g. leakage). In this paper we 
point out that such a competitive effect is essential e.g. 
for the solid-state qubit implementation using excitonic 
(charge) states in quantum dots (QDs) with com- 



putational states defined by the absence (|0)) or pres- 
ence (|1)) of one exciton in the ground state of the dot, 
operated by resonant coupling to laser light. In such a 
system, Markovian decoherence (exponential damping) is 
related to exciton recombination on 1 ns timescale, while 
strongly non-Markovian effects result from lattice iner- 
tia. The interplay of these two decoherence mechanisms, 
favoring opposite strategies (fast vs. slow gates), leads 
to a kind of trade-off, resulting in optimal speed of gate 
for most reliable operation on the qubit. 

To investigate the effect, we analyze decoherence in 
a general spin-boson system P, jl^], and find decoher- 
ence to be strongly non-Markovian, fitting into "minimal 
decoherence model" . Assuming additional, Markovian 
damping, we obtain the trade-off formula for the error 
caused by decoherence, averaged over input states of the 
qubit 
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(actually, in the solid state example, singularity of the 
first term is lifted by the presence of upper cut-off). The 
constants 7m, nM express the strength of the Markovian 
and non-Markovian decoherence, respectively. We deter- 
mine their values for the exciton in a QD with typical 
parameters and show that the gate duration leading to 
minimal overall error is of the order of 1 ps. Since level 
spacing in such a system allows even 100 fs gating |^ the 
restrictions imposed by non-Markovian phonon mecha- 
nisms are decisive for the gate speed. Our result has 
both practical as well as general implications for QIP: 
(i) it suggests the proper direction in research towards 
semiconductor implementation of QIP (ii) it provides a 
non-trivial dynamical strategy to minimize decoherence 
in quantum systems, for a class of reservoirs. 

Minimal decoherence in the spin-boson model. We con- 
sider a qubit described by means of the spin-boson model. 
The Hamiltonian H" of the system plus reservoir is 



Hi 



Hs{t) + HsB. + Hr, 



(2) 
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where S is our qubit system, R is (bosonic) reservoir. 
Hg — ( 1/2) ficT J, is the self Hamihonian of the qubit, 
Hs{t) = (l/2)e(t)(e*"*cr+ + e-^^^V-) is the gate Hamil- 
tonian (in rotating wave approximation), e(t) being the 
shape of the pulse. Hji = J^'^kO-kil. is the self Hamil- 
tonian of the reservoir, where Wfc is the boson energy. 
Finally, Hsb, is the interaction with reservoir 



HsR — 



\Y1 fk^k 



H.c. 



(3) 



where is the coupling constant for the mode fc. 

In the interaction picture with respect to Hg we get 
H' = e{t)ax + HsR + Hr. Let us now represent H' in the 
basis of eigenstates of the total system (dressed states). 
This is achieved by the unitary operation U = |0)(0| ® 
W - |1)(1| ® W\ where W = eMY.Ul/^k)ak - H.c] 
and |0), |1) are the eigenvectors of ctz- In this basis the 
Hamiltonian H = U^H'U up to a constant has the form 



H = 



^e{t)ar cos ( 2i — aj, + H.c. ) 
2 \ k '^'^ J 

+ ^eit)ay sin (^i ^ + H.c. 



(4) 
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To the first order in fk/^k we obtain 
H ~ e(i)cr^ + e{t)ay ® [ iY^ —a) + H.^ 



Hr. (5) 



where []av denotes averaging over the states ■0 and 
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A complete minimization of 5 would require full opti- 
mization of the pulse shape. However, in order to demon- 
strate the idea of the trade-off in simple terms, let us 
restrict the discussion to qubit rotations performend by 

Gaussian pulses, e{t) — a/(V27rTg)e~^ . Here Tg is 

the gate duration, while a is the angle determining the 
gate, e.g. a = -I is VNOT, while a = tt \s (Jx (bit flip). 
The functions |F±(ti;)p may be written as 



(9) 



Although the original minimal decoherence model with 
its ~ us^ dependence appears in many physical situations, 
other characteristics are obviously possible. In general, 
as may be seen from (jSj), JT)) and for a spectral den- 
sity Rioj) ~ oj" the error scales with the gate duration as 
T~"+^ and Tg""^^ at low and high temperatures, respec- 
tively. Therefore, for n > 2 (typical e.g. for various types 
of phonon reservoirs •ill) the error grows for faster gates. 
Assuming the spectral density of the form Riuj) — Rquj^ 
(see the QD example below), we obtain 



We define the qubit in terms of the dressed states [I^ ■ 
With such a choice, the reservoir is decoupled and the 
interaction term is present only during gate operation. 
Thus our system may be reduced to the minimal deco- 
herence model: decoherence is not present at all, if the 
qubit is not active. Now we want to calculate the result- 
ing error and relate it to the speed of gate. Since the 
regime is strongly non-Markovian we solve the Master 
equation in the Born approximation and compute the fi- 
delity F — 1 — 6 of a. single-qubit operation (see Q for 
details) . 

The error, averaged over the initial qubit states can be 
represented as the overlap of two functions 
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duj 



R{lo)S{uj). 



(6) 



Here R{u!) is the spectral density of the reservoir 
R{uj) = [uBiio) + 1] ^ [d{cOk -io) + <5(c^k + u;)] \fk\^ 

k 



(7) 

where nB(w) is the Bose-Einstein distribution. The func- 
tion S{uj) fully represents the spectral characteristics of 
the system and is given by 

s{lj) = [(V'|yrt|^)_(^|yt|^)(^|y|^)]^^ (8) 

« l(|F_(.;)P + |F+(c.)p), 



at T = 



(10) 



This leading order formula holds for S ^ 1. Also, if we 
introduce the upper cut-off, the error will be finite even 
for an infinitely fast gate [ijl (see Fig. 

Trade- ojf between two types of decoherence. As we have 
shown, in our model the error grows as the speed of gate 
increases. This could result in obtaining arbitrarily low 
error by choosing suitably low speed of gates. However, 
if the system is also subject to other types of noise this 
becomes impossible. Indeed, assuming an additional con- 
tribution growing with rate 7m, the total error per gate 
is 



e 7nM , 

— 5- + 7MTg, 7nM 



1 2 1 , ^ 

-a Rq, 7m = — , (11) 

3 Ti- 



where Ti- is the characteristic time of Markovian decoher- 
ence (recombination time in the excitonic case). As a 
result, the overall error is unavoidable and optimization 
is needed. The above formulas lead to the optimal values 
of the form (for T 0) 



i n — 



3 /2a2i?o^ 



1/3 

for Tg ^ I ^a'RoT, ] . (12) 
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We will see that the trade-off is present also for finite 
temperatures. 
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Example: exciton in a quantum dot (QD). Let us now 
estimate the error magnitude for the specific semicon- 
ductor QD qubit implementation 0. The reservoir is 
then constituted by phonons, the most important branch 
for our present purpose being the longitudinal acoustical 
(LA) one, characterized by linear dispersion, ojk — ck, 
coupled via deformation potential (DP) JJj. In fact, the 
relatively simple model is very accurate for a descrip- 
tion of the QD system for ~ 1 ps timescales relevant here, 
as confirmed by the excellent agreement between the the- 
oretical calculations and experimental results 
This is due to the fact that neither the high-frequency 
optical phonons nor direct or phonon-induced leakage 
to higher states contribute considerably to the decoher- 
ence. Also effects related to piezoelectric coupling to LA 
and TA phonons may be neglected in weakly piezoelec- 
tric systems (e.g. GaAs). Moreover, the qubit control 
is all-optical, eliminating the need for additional noise- 
inducing device structures. On the other hand, details of 
the QD structure (shape, stress, composition) may lead 
only to quantitative modifications of secondary impor- 
tance. 

The non-Markovian error has a dynamical origin and is 
related to a which path trace left by exciting the phonon 
modes rather than to an influence of noise. The spec- 
tral density is uniquely defined by the lattice response 
characteristics: the mode-dependent coupling strength 
and the density of states. Due to fundamental restric- 
tions (global charge neutrality, translational invariance), 
the frequency dependence is always strongly super-ohmic 
|l5j . On the other hand, the approximate momentum 
conservation holding for weakly confined carrier states 
leads to exponential cut-off of carrier-phonon interaction 
at high frequencies [31 ■ In the s pec ific case of DP cou- 
pling, the coupling constants are |12| 




(13) 



where v is the volume of the unit cell, p is the crystal 
density, (Te,h are deformation potential constants for the 
electron and the hole (material parameters are taken as 
in ^^) and ^j^'^^ are formfactors for the correspond- 
ing wavefunctions (approximated by Gaussians, on the 
grounds of numerical diagonalization in parabolic con- 
finement 
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FIG. 1: The total error for a. a = ti/2 rotation on a QD qubit, 
for T = (solid lines) and T = 10 K (dashed lines), for two 
dot sizes {Ih = O.SZe, Iz ~ 0.2le). The Markovian decoherence 
times are inferred from the experimental data [T3 . Inset: 
Spectral density of the phonon reservoir R{iL>) at these two 
temperatures and the gate profile ^(tj) for a — n/2. 



the QD system, may be interpreted as "dressing" the 
localized carriers with coherent lattice deformation field 
(cf. ^Td\). This kind of dynamics has been described 
in the limit of very rapid gating pulses [Til ITsI l and has 
been shown to lower the degree of coherent control over 
the system and to destroy the coherence of polarization 
oscillations 0, 0| . Due to relatively large dot sizes and 
low sound speed the cut-off frequency (~ 1 meV) is lower 
than the transition energy between different carrier states 
in a small, self- assembled dot and the corresponding 1 ps 
times are experimentally accessible. 

For large enough gate duration only the low-frequency 
sector contributes and the coupling H13|l may be approx- 
imated by 



-(o-o - fJh) 



leading, according to Q), to 



R{lu) ~ i?o^^^ Rq = 



2gvc 



(15) 



(16) 



Hence, Eqs. ((T^ are applicable and one finds for the 
specific material parameters of GaAs, 



■J I 



(e,h) 



(14) 



Tg = a^l^XAl ps, (5min = a^/^0.0035. 



(17) 



where Zc.h and Iz are the wavefunction widths in the dot 
plane, for electron and hole, and in the growth direction, 
respectively, and fcx and kz are the corresponding com- 
ponents of the phonon wavevector. 

Note that the most effectively coupled phonon modes 
correspond to wavelengths around the dot size. For gap- 
less bosons with linear dispersion, their frequency de- 
termines the reservoir memory times and sets up the 
timescale of non-Markovian effects which, in the case of 



The full solution within the proposed model, taking 
into account the phonon cut-off for an anisotropic shape 
ij-z < h,h), according to Eqs. ifHI) . ((T^ . and allowing 
finite temperatures by numerically calculating the spec- 
tral density Q, is shown in Fig. The size-dependent 
cut-off is reflected by a shift of the optimal parameters 
for the two dot sizes: larger dots admit slightly faster 
gates and lead to lower error. Interestingly, the trade-off 
becomes more apparent at nonzero temperature. 
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The presence of the upper cut-ofF could suggest apply- 
ing the dynamical decoupling (DD) technique 0| to di- 
minish decoherence. However, for high frequencies many 
other mechanisms of decoherence become relevant. On 
the other hand, combining the bounded-control version 
of DD [l^ with the optimization proposed here might 
lead to some reduction of the resulting error. Such tech- 
niques might also be useful for eliminating other sources 
of noise, not included in our discussion. 

Another possibility to reduce the decoherence effect in 
the QD system might be to encode qubits in decoherence- 
free subspaces jl9|. However, the phonon wavelength cor- 
responding to the optimal range of gating times is com- 
parable to the single dot size, precluding the necessary 
collective interaction with the whole QD system. Thus, 
for decoherence effects which do not involve real transi- 
tions and taking into account the feasible system geom- 
etry and actual nature of phonon coupling, only a small 
decrease of the minimal decoherence may be expected. 

In conclusion, we have exhibited the trade-off between 
two opposite types of decoherence: usual Markovian 
damping and dynamically induced non-Markovian deco- 
herence for a realistic super-ohmic reservoir. To protect 
the qubit, opposite elementary strategies are needed: fast 
and slow gates, respectively. The minimization of the 
overall error leads to optimal speed of gate. We have 
shown that the trade-off is present in a semiconductor 
implementation of quantum information processing. The 
Markovian error is caused by recombination, while the 



non-Markovian one occurs if the gate operation is not 
adiabatic with respect to lattice modes. We have evalu- 
ated the minimal error in this case for a single qubit gate, 
showing that the two processes indeed compete. The op- 
timal gating time 1 ps) sets up the limit beyond which 
any further gate speed-up is unfavorable. Even at this 
optimal point the trade-off gives rise to a significant er- 
ror. For two qubit gates involving single-qubit rotations 
(as proposed e.g. in j^), the present result gives a rough 
lower bound for the error, which is of 1-2 orders of magni- 
tude higher than that admitted by fault-tolerant schemes 
Q known so far (~ 10~^). However, possible improve- 
ments of the latter schemes cannot be excluded. It follows 
also that diminishing the responsible constants 7m (e.g. 
by elimination of radiative losses |^) and (by op- 
timizing the system parameters) is the most important 
task towards semiconductor implementation of a quan- 
tum computer. It is also important to explore whether 
the same effect can occur in other implementations, as 
well as to what extent the error avoiding techniques may 
be helpful here. 
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